Computing Borel's Regulator I 



Zacky Choo, Wajid Mannan, Ruben J. Sanchez-Garcia, Victor P. Snaith 



September 4, 2009 



Abstract 

We expand Hamida's integral for the Borel regulator in odd dimension as an infinite 
series and give a computer algorithm for dimension 3. 



1 Introduction 

Let Cf{s) denote the Dedekind zeta function of a number field F. The analytic class 
number formula ([i9_ p. 21) gives the residue at s = 1 in terms of the order of the 
class-group of Of, the algebraic integers of F, and the Dirichlet regulator Ro{F). Let 
dp denote the discriminant of F. In terms of algebraic K-groups of Of the class-group 
is equal to the torsion subgroup Torsive (Of) of Ko{Of) and the formula has the form 

^ , , 2'-'+'-^7T'-^Ro{F)\TorsKo{OF)\ 
ress=i Cf(s) = ^ 



|Tors Ki(OF)\Vd^ 

The Dirichlet regulator Ro{F), which is a real number, is the covolume of the lattice 
given by the image of the Dirichlet regulator homomorphism f [19] p. 25) 



R°F --Of =Ki{Of) 



jri+r2 — 1 



Here ri and 2r2 denote the number of real or complex embeddings of F respectively. 
Equivalently, by Hecke's functional equation ([12], [19] p. 18), Cf{s) has a zero of order 
-I- r2 — 1 at s = 0. Let ^|-(so) denote the first non-zero coefficient in the Taylor series 
for (^F Sit s — So- Therefore at s = the functional equation yields 

a(0) = lim ^^^'^ - MF)\^orsKoiOF)\ 



s_o s'-i+'-2-i |TorsJi'i(OF)| 

This form of the analytic class number formula prompted Lichtenbaum [ll] to ask: 
Which number fields F satisfy the analogous equation for higher-dimensional algebraic 
K-groups 

.,,.^,^.MF)\K^2r{OF)i 

' |Tors7fi_2.(OF)| 
for r = —1, —2, —3, . . . and some integer e? Here Rr{F) is the covolume of the Borel 
regulator homomorphism defined on Ki-2riOF) and to which we shall return in §2. 
This identity has become known as the Lichtenbaum conjecture and is known to be 
true in many cases; for example, there are several proofs for abelian extensions of the 
rationals ([H], [IS], [16], [l^). 

In the special case of the integers the conjecture takes the form 



RriQ) |TorsA'i_2.(Z)l 

for r = —1, —2, —3, . . . and some integer e. Here (^*{r) is the leading coefficient in the 
Taylor series for the Riemann zeta function ((z) aX z = r. As mentioned above, this 
and the analogous relation in the cyclotomic case, is known to be true. Moreover when 
r = —1 the value of e is known in the cyclotomic case by the main result of Ion Rada's 
thesis [14]. In addition the value of |Tors /('i_2r(C'F)| (for example, see [16], [17]). The 
nett result is that, in order to compute |_ft'2(C'F)| in the cyclotomic case it suffices to 
have merely approximate values for Cf(^1) ^.nd of the Borel regulator R-i{F). 
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This discussion provides the motivation for this paper. For example [TS] the Kummer- 
Vandiver conjecture, which originated in an 1849 letter from Kummer to Kronecker, 
would be implied by the vanishing of |_K'_2r(Z)| for r = —2, —4, . . .. Although we are 
only concerned here to design an algorithm to calculate an approximation to the Borel 
regulator in dimension three in the long run we hope that our ideas may be capable of 
implementation in higher dimensions to calculate orders of even dimensional K-groups 
of rings of algebraic integers. 

2 The Borel regulator map 

Let _F be a number field and Of its ring of integers. F have [F: Q] different embeddings 
(Tj into C. Some of them, say ai, . . . , , are real and the rest come in conjugate pairs 
(Jri+i, . . . , o"ri+r2 ) o'n+i, . ■ . ,o'ri+r2- The Borel regulator maps are homomorphisms 

K2n+i{OF) ^M."" (n>l) (1) 

from the odd algebraic A'-theory groups of Of to R'*" , where d„ is ri + r2 if n is odd 
and r2 if n is even. Given an embedding a: F ^ C, write ct* for the induced map in 
K-theory. Borel defined homomorphisms 

B„: A'2n+i(C) (2) 

such that the map Q is the composition 

K2u+i{Of) ^ K2n+1 (C) ^ e C . 

The map Bn is in turn defined as the composition of the following homomorphisms. 
The algebraic K-iheoiy groups can be defined as the homotopy groups 

A'2n + l(C) =7r2n + l(BGL(C) + ), 

and hence there is a Hurewicz map into homology with integral coefficients 

K2„+l{C) H2n + liBGL{C) + ) ^ J/2n + l(GL(C)) . 

To obtain the isomorphism above, recall that the plus construction does not change 
the homology and the homology of a discrete group G is the homology of its classifying 
space EG. Via Suslin's stability results [S] , if A'' > 2n + 1 then 

i?2„+l(GL(C)) ^ H2n+liGLN{C}) . 

Therefore, it suffices to define a map 

6„:H2n+i(GLiv(C))— (3) 

and precompose with the homomorphisms above. The map ^ is called the universal 
Borel class and it is actually given by a cocycle in the continuous cohomology group of 
GLiv(C). The image of 6„ lands in (27ri)"~^R which we identify with R. The original 
construction, due to Borel and using continuous cohomology, is described in [3]. Borel 
also proved that the image of the map ([1} is a lattice in R"*" , whose (co)volume is called 
the Borel regulator, written Rn{F), as mentioned in §1. 

3 Hamida's integral 

By a result of Hamida jB] , the universal Borel class has the following description as an 
integral of differential forms. Fix n > 1 and Xq, . . . , Xn G GLn{C) with A > 2n + 1 
(ie. in stable range). Let A„ be the standard n-simplex 

A„ = |(a;o,...,x„) G R"+' | > 0, ^ i« = ij. 

Define for every point x = (xq, . . . , x„) G A„ 

u(x.) = xoXqXo + . . . + x„X*X„ , (4) 
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where X* denotes the conjugate transpose of X. Thus is a matrix of 0-forrns (com- 
plex functions) on the n- manifold A„, and vijx) is positive definite hermitian and in 
particular invertible. Consider the matrix of differential n-forms {u~^du)^^^ and define 

(^(Xo,Xi,...,X„) = Tr / {p-^dv)"'' (5) 



the trace of a matrix of integrals of differential forms. 
Theorem 1 (Hamida [S]) Letn>0. The map 

(Xo,...,X2„+i) - ^=l):^^iXo,X„...,X,^+,) (6) 

induces a cocycle m the continuous cohomology o/GLjv(C) which is the universal Borel 
class b„. 

Remark 2 The cocycle (0) is homogeneous and unitarily normalised, that is, 

ip{Xog,...,Xr,g) = <^(Xo,...,X„) for all g G GLiv(C), (7) 
ifiiuoXo, ■ ■ ■ ,u„X„) = ip{Xo, . . . ,Xn) for all G l7jv(C). (8) 

In particular, we can assume Xn ~ Id by Q), and all the Xi to be positive definite 
hermitian matrices by via the polar decomposition: every invertible matrix X can 
be written as X = UP where U is unitary and P is positive definite hermitian. Indeed, 
X* X — P* P; compare with 

4 The Infinite Series Formula 

Our goal is to make the computation of Hamida's integral (O explicit. Namely, we will 
transform the integral into an infinite series which we can arbitrarily approximate. 

In order to integrate an n-form on an n-dimensional manifold (with boundary in 
this case) we need to express the integrand in terms of n cordinates rather than n + 1. 
We have a homeomorphism from the n-simplex given in R" by 

n 

A" = {(yi, 2/2, . . . , i/n) i ?/s > for aU i , ^ y^ < 1} 
to the n-simplex in R""*"^ given by the map 

n 

iyl,■■■,yr^) (i^X] . . . ,yn). 

Therefore, in terms of the j/-coordinates, we have 

u = X^Xq + yiiX^Xi - Xo*Xo) + . . . + yn{X:X„ - XSXo) 

~ X*{{X*) ^XqXoX^^ +yi{{X*) ^XiXiXn^ — {X*) ^X*)XoX„^) 

+ .. . +y„(l - {X-*)-^XSXoX-^))X„ 

and 

du = x*{dyi{{x:y'x:x^x-' - {X*)-'XSXoX-') 

+ ...+ dyn{l - {X:^)-'X5XoX-'))X„ . 
Therefore v = X^v'Xn, dv — X^di/'Xri and 

Tr / (p^^du)" = Tr / {{u')-^ du')" . 

Next, we perform a change of variable by means of the map 

r : [0, 1] X A""' > A" 

given by 

T(t, Sl, S2, . . . , Sn-l) — {sit, S2t, . . . , Sn-lt, 1 — t). 
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Since sit + S2t + . . . + Sn-it+l — t — 1 — — si + S2 + ... + s„_i) lies between zero and 
one, the image of T lies in the rz-simplex. Furthermore for each non-zero value of t the 
corresponding horizontal (n — l)-simplex is mapped homeomorphically onto its image 
while {0} X A"~^ is collapsed to the vertex (0, 0, ... , 0, 1). The Jacobian of T is equal 
to Therefore when f / we have J(T) > if n is even and J{T) < if n 

is odd. For any compact n-manifold M with boundary in [0, 1] x A"^^ having image 
T{M) in the n-simplex we have (g] p. 28) 

J^j if ■ T)(t, si, . . . , s^.i){-l)"r-^dtdsids2 . . . ds^-i 

= It{m) fiyi:---^y^)dyi---dyn. 

We shall be interested in the integral Trace {{u')^^dv')" which may be written 
as the limit of integrals over manifolds of the form T{M) as they tend towards the 
n-simplex. Therefore we can compute this integral as a limit of corresponding integrals 
over M C [0, 1] x A"~^, provided that this limit exists. However the integral over M, 
involving the Jacobian of T, is merely the integral over M where • T is written in 
terms of t,si, . . . , Sn-i and di/ is also computed in these coordinates. Thus we have 

v' ■ T = (X*) ^X^XoXn^ +tsi{{X*) ^XiXiX„^ — {X*) ^X^XoX:^^) 

+ . . . + (1 - t){l - {X*)-'Xt,XoX-') 
= {X^)~'^XoXoX^^ + tsi{{X*)^^XiXiX^-^ — (X*)"^XoXoX^^) 

+ . . . + 1 — t — {X*) ^XoX()X„ ^ + t{X*) ^Xt)XoX„ ^ 
= tsi{{X^) ^ XiXiX„ ^ — (X^) ^ Xq XoX„ 

+ ... + l-t + t{X:i)-'XSXoX-' 

= 1 + tA{si,S2, . . . , Sn-l) 

where 

A{S1,S2, ■ ■ ■ , Sn-l) 

n-1 

= {{X:)-^XoXoX-' + Sj{iX*n)-'x;X,X-' - (X*n)-^X*oXoXn'). (9) 

Observe that, if we assume that the Xi's are close together, then for all points of A"~^ 
the matrix A{si, S2, ■ ■ ■ , Sn-i) is close to zero. Under this assumption we may express 
the integral as a convergent infinite series in the following manner. 
As M varies, we must examine the integral 

Tr / {{l + tA)~\dtA + tdA))" 

J M 

where dA = J2]Zi dsjUXiy X* XjX'^ - {X'nY^ X^XoXn^). Hence we obtain 

Tr/^,^ {{l + tA)-\dtA + tdA)Y 

= Tr/^^ (1 + tA)~^dtA{(l + tAyhdA)'^-^ 

+ Tr /^^ (1 + tA)-hdA{l + tA)-^dtA{{l + tAyhdA^-^ 
+ Trjj^j ((1 + tA)-^tdAf{l + tAy^dtA{{l + tAyhdA)"-^ 



+ TtJj^j {{1 + tA)-hdA)"-^(l + tA^dtA. 
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since the integral is unchanged by cyclic permutation of factors and because A{s) and 
1 + tA{s_) commute. Next we observe that if Wi, . . . , Wn are m x m matrix- valued 
functions then 

TTace{WidxiW2dx2 . . . W„dx„) = (-l)""'Trace(Ty2rfs2 . . . W„dx„Widxi) 

because the trace of a product of n matrices is invariant under cyclic permutations but 
the n-form dxi . . . dxn changes by the sign of the n-cycle, which is (—1)"^^. Accordingly 
each integral in the formula for 

Trace [ {{1 + tA)-'^{dtA + tdA))" 

J M 

is equal to ( — 1)"^^ times the one below it. Together with the barycentric subdivision 
results of the next section this observation implies that 

V9(Xo,Xi,...,X2„) =Tr / {u-^dvf' ^Q. 

For the rest of this section we shall assume that n is odd. In this case 
Trace /[o,i]xA"-i ((1 + tA)-\dtA + tdA)r 
= n Trace /[oi]xA"-i + tAy^dtA^il + tAy^tdAY'^ 
= n Trace 1 r-^dtA{l + tA)-^{{l+tA)-^dAy'~^ 

= n Trace ^j^^,.-! f'-^dtAil + tA)-'^dA{l + tAy^dA.. . (1 + tA)-^dA. 

Also we observe that the assumption that the XiS are close together means that 
through the domain of integration the matrix A is very close to zero. Recall the 
geometric series formula for a matrix A [71 5.6.16]: if || • || is a matrix norm and \\A\\ < 1 
then I — A is invertible and YlT^o"^^ ~ ~ ^)~^- matrix norm on M]v(C) is a 
vector norm which satisfies ||^-B|| < ||A||||_B||.) Therefore we may expand the integral 
as a series, taking care not to commute matrix factors which may not commute, to 
obtain the following formula: 

Trace J^^^^^^^„_, {{1 + tA)-' {dtA + tdA)^ 

= n Trace /j^ f^'^dtAil + tA)-^dA{l + tA)-'^dA . . . (1 + tA)-'^dA 

= nTrace /[o.i, ^ a-i Eo<™. + 1)^"' dtA{tAr^ dA . . . {tAr^-^dA 

= nTrace /,g ^^^,,-1 E„.>o (-l)-''"i'--™"-i'"'mit2i(™i'--'"'.-i)+"-idfyl'"idyl . . . 

where rn{mT_, . . . , m„_i) = mi + . . . + m„_i. 

In §5 we shall develop a barycentric subdivision algorithm to reduce to the case 
where the Xi's are close together so that < 1 and then we shall be able to evaluate 
the infinite series expression 

nTr / V (-ip(™i--™"-i'-^mit-''"i'--™"-i)+"-'dtA™idyl...A'""-idA, 

i[0,l]xA"-l 

(10) 

where m{rn\, . . . , m„_i) = m\ + . . . + m„-i. 

4.1 Case n = 1 

In this example the map T map be taken as the identity map and M — [0, 1]. We 
obtain 

TrJ^^iu^^du) ^TrJ^{l + tA)-^Adt 

= TrE„>o/o(-l)"i"^"^rfi 
= Tr(log(l + A)). 
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4.2 Case n = 3 

For n — 3 we can simplify (|10p further. We have 
Trace J^^^ {{1 + tA)"'^ {dtA + tdA)f 

= 3Tia.ceJj^jJ2m,,m2>o {-l)""''^"'^~''mit"'^+"'^+^dtA"'^dAA"'^dA 

and because t^i+^^+a^^ jg ^ ^^^^^^ 1-form the effect of switching A'^^dA with A'^^dA 
is simply a change of sign 

Trace /j^^ ((1 + + td^))^ 

= -3Trace f^^T, >n (-l)'"i+'"'"'mit™i+'"2+2df^'"2dAA'"idA 

J A/ 'mi ,7712 >0 ^ ' 

= -STrace T^Y] (-l)"'^+"'^-'^m2t"'^+"'^+^dtA"'^dAA""^dA. 
This yields a symmetrical formula 

Trace /^^^^ {{1 + tAy^ {dtA + tdA)f 

= |Trace/^^5:;^^_,^^>o (mi - m2)t'"i+™2+^dfyl"idylyl'"2dyl. 

Now let us integrate the formula 

Trace /^^ {{1 + tAy^ (dtA + tdA)f 

= |Trace/^^Emi,™,>o (-l)'"^+'"^"'(mi - m2)t'"i+'"^+'dM'"idAA'"^dA 
= 3Trace/^^Emi>m2>o (-l)'"^+'"^"'(mi - m2)f'"i+'"^+'dtA'"idAA'"2dA 
^STraceLaV ^ (_l)'"i+'"2-ii!2izi:!jI^mi^^^m2^^ 
Hence the intregral for n — 3 and ||y4|| < 1 is 



3Tr 



Ja^ ^ mi + ma +3 ^ ^ 

mi >m2 >U 



5 Barycentric subdivision 

In order to use the infinite series formula in Section 14] we need \\A\\ < 1 for a matrix 
norm |j • ||. The idea is to apply the iterated barycentric subdivision operation to the 
simplex represented by {Xq, ■ ■ ■ , Xn) until the associated matrix A of each subsimplex 

satisfies \\A\\ < 1. 

5.1 Definitions and main properties 

Suppose that X is a normed real vector space. Define the simplex spanned by xq, xi, . . . , Xn G 
X as the subset 

A[a;o, ...,x„]^ {toxo + ...+ t„x„ \ U G R, U > 0, ^ = 1} . 

i 

The barycenter of A[2;o, . . . , Xn] is the point 

b = — — r(^o + • ■ ■ + x„) . 
n + 1 

The barycentric subdivision of A[a;o, . . . , a;n] (cf. [18]) is the following decomposition 
into (n + 1)! n-simplices. For each permutation a G En+i we have the n-simplex 

Sct = AfXo-jo), a;o.{0)o-(l), a^(T{0)<T(l)CT(2), • • • , ^cy{0)cy{l)...a{n)\ , 

where we write 

m 

the barycenter of A[xij^, . . . ,Xi^]. We define the sign of s^ as sgn((T), the signature of 
the permutation. 
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Remark 3 For these definitions to make sense, we only need X to be a convex cone 
in a real vector space. For example, the set of all positive definite hermitian matrices. 

The barycentric subdivision reduces the diameter by a factor of 

Proposition 4 Suppose that A ts an n-simplex in a normed real vector space spanned 
by xo, ■ ■ ■ ,x„. Consider its diameter diam(A) = sup{||a; — y\\ : x,y £ A}. Then 

(1) diam(A) = maxij \\xi — Xj\\; 

(2) if A' is any n-simplex in the barycentric subdivision of A then 

diam(A') < diam(A) . 

^ ' - n+1 ^ ' 

Proof: Given a,b £ A, suppose b — tiXi as above. Then 

< max \\a ~ Xi\\ . 

i 

A similar argument shows \\a — Xi\\ < maxj \\xi — Xj\\ and the first claim follows. For 
(2) see, for instance, |18l pp. 124-5]. 

Remark 5 The barycentric subdivision procedure can be iterated, by applying it to 
each of the resulting n-simplices. Hence we can talk of the kth barycentric subdivision 
of a simplex, for each k > 0. 

The second key property is that barycentric subdivision can be applied and induces 
the identity in group cohomology. Let G be a group, R a commutative ring and M 
an TiG-module. Consider {F„,dn) the standard resolution of RG Each F„ is the 
free /JG-module generated by (n + l)-tuples {go, . . . , gn) with G-action g{go, ■ ■ ■ , gn) = 
(55o,---,Sffn) and 

n 

du{go, ■ ■ ■,g,i) = ^(-l)'(5'o, . . . ,5,, . . ■,g,i)- 

If we define C„ = F„^rgM , then the homology groups HniG; M) are the homology 
groups of (G*, 9* ® 1m). On the other hand, if we let G" = YioxanGiFn, A4), then the 
cohomology groups H"{G;M) are the cohomology groups of [C* ,5*) where 

n+l 

'5"(")(5o, • ■ ■,gn+i) = ^(-1)'q(3o, ■ ,fln+i) for aU a G G". 

The barycentric subdivision is a chain map S: F„ —> Fn defined on generators as 
S(go, . . . ,5n) = ^ sgn(a)(g^(o), . . . ,p^(o),...,a{n+i)), 

where 

= — (ffii + . . . +gi,J G RG. 
m 

(From now on we assume char(_R) = so that 1/m £ R for every m > 1.) 

5* extends to (co)chain maps 5": G„ ^ C„ and S: G" G" (we use the same letter 
S). For example, 

S{a){go, . . . ,gn) ^ ^ sgn((j)a(5^(o), 51^(0),.. .,<T(n+i)) for aU a G G". 

Remark 6 The barycentric subdivision can be iterated, simply by = S o . . . o S , 
k > 1 times. 

It is standard and not difficult to show that these chain maps induce the identity in 
(co)homology. In particular, if a G G" is a cocycle, then so is S{a) and the represent the 
same element in H"{G; M). Recall that Hamida's function is a cocyle representing 
(up to a coefficient) the universal Borel class in H"{GLn{C);C). Hence S{ip) also 
represents the universal Borel class. The same holds if we further subdivide. 

Proposition 7 Let k > 1. The cocycles ip, as defined in (3)), and S''{(p) represent the 
same element in iif"(,„i(GI/jv(C); C). 
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5.2 On matrix norms 

Let II • II a matrix norm, that is, a vector norm in Mjv(C) which satisfies ||j4_B|| < 
II j4|| ||_B|| . Our main example wiU be the spectral norm [7] 5.6.6] 

||A||2 — max : A is an eigenvalue of . 

This norm satisfies ||A*||2 = ||^||2 and ||[/v41/||2 = ||A||2 if U and V are unitary. If 
A is hermitian then ||A||2 = max{|A| : A is an eigenvalue of A}. In this case, all the 
eigenvalues are real and we write them in increasing order as 

Xrmn{A) = Xl{A) < A2 (A) < • • ■ < A„(A) = A^^^A). 

If A is positive definite hermitian (if and only if \min{A) > 0) then ||yl||2 ~ Xmax{A) 
and P~^||2 = l/X,nin{A). 

Theorem 8 (Weyl) Let A and B be hermitian matrices of dimension n. For each 
k = 1, . . . ,n we have 

Xk{A) + Ai(B) < Afc(A + B)< Xk{A) + A„(B). 

For a proof, see 4.3.1 in [7]. 

In particular, if Ai, . . . , A„ are hermitian matrices then 

Xmin{Ai-\ + A„) > Xmi„{Ai) -\ + A™,i„(A„) and (12) 

Xmax {Al + ■ ■ ■ + An) < Xmax{Al) + ■ ■ ■ + Xmax{An). (13) 

5.3 Upper bounds 

Let Xo, . . . , Xn G GLn{C) and Si, . . . , s„-i G R with Si > and Si < 1. Define, as 
in §4, 

n-l 

A = {{X:)-'XoXoX-^ ^i) + Y^s, {{X:)~^X;X,X-^ - iX:)~'XoXoXn') . (14) 

3=1 

Lemma 9 Let \\ ■ \\ be any matrix norm such that \\A*\\ = ||j4|| for all A G M„(C). 
Then, for any matrices A,B,C G M„(C), we have 

(a) WA'A - B'BW < ||A - B|| ■ (P|| + ||B||), 

(b) \\C*A*AC - C*B*BC|| < ||C|p ■ \\A - B\\ ■ {\\A\\ + \\B\\), 

(c) \\C*A*AC - 111 <\\C\\^-\\A- C\\ ■ {\\A\\ + \\C\\). 
Proof: 

(a) \\A*A - B-BW = ||(yl* - + B*{A - B)|| < p - B\\ ■ \\A\\ + \\A ~ B\\ ■ \\B\\ 

< P-B||-(||A|| + ||B||). 

(b) \\C*A'-AC - C*B*BC|| = ||C* {A* A - B*B) C\\ < \\Cf ■ WA'A - B*B|| and use 
part (a). 

Part (c) follows from (b) when B = C^^ . 
Using Lemma [9] we obtain 

(n — l 
\\X„ - Xoll • (||X„|| + IIXoll) + ^s,||X, - Xoll ■ (||X,|| + ll^oll) 

^ n — 1 

< ||X-^f ■ max ||X, -Xoll ■ ( ||Xn|| + ||Xo|| (||X,|| + ||Xo||) 

< ||X-if ■ max ||X, -Xoll ■ M|X„|| + max ||Xj || + 2||Xo | 

l<i<n \ l<j<Ti — 1 

Proposition 10 For any matrix A as in ji^P ^'^'^ "-^V matrix norm satisfying \\X* \\ = 
||X|| for all X G M„(C), we have 

Pll < ll^n'll'- max ||X, -Xoll- f||X„||+ max || JV, || + 2||Xo || V (15) 

l<i<n \ l<j<n—l J 
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Corollary 11 Let d = diam(A[Xo, . . . ,X„]), M = maxi ||Xi|| and m = maxi \\X^ ^\\. 
Then 

\\A\\ < AdMrn. 

If the matrices Xi axe positive definite hermitian, we can consider iterated barycen- 
tric subdivisions of A[Xo, . . . , X„] (cf. Remarks |3] and [2]) . Then the bound above, with 
respect to the spectral norm, decreases arbitrarily. 

Proposition 12 Suppose that Xa, . . . , Xn are positive definite hermitian matrices span- 
ning a n-simplex A — A[Xo, . . . ,X„]. Let \\ ■ \\ be the spectral norm and define 

d = diam(A) = max II Jfi — Jfj II , 

M = ma.x{\\Xi\\} = ma.x{Xmax{Xi)} , 

i i 

m = max{||X,"^||} = max{l/Amin(^i)} ■ 

i i 

Let A[Xo,...,X„] he any n-simplex of the k-th iterated barycentric subdivision of 
A[Xq, . . . , X,i]. For si, . . . , s,i > and "^Si < 1, define 

n-l 

A = (^{Ky^XoXoX-' - /) + ^ ((X 

) ^j^n ^ i^n) XoXoXn 

J = l 

Then 



k 



||I||<(-^J .4rfMn^^ 

This guarantees that, for k large enough, all the matrices A as above will satisfy 
||j4|| < 1 (the bounds d, M and m only depend on the initial matrices Xo, . . . ,Xn). 

To prove the proposition, we need first a bound for ||X|| and ||X~^|| for an arbitrary 
point in A[Xi, . . . ,Xn\. 

Lemma 13 Let Xq, . . . , X„ be positive definite hermitian matrices and X G A[Xo, . . . , Xn] 
Suppose that X = X]r=o ^i-^i with ti > 0, X^f^o ~ ^- Then 

\\X\\ < ^t,||X,|| and 

i=0 

n 

In particular, \\X\\ < max^ ||Xi|| and ||X^^|| < max^ H^j^^H. 

Proof of Lemma: The first inequality holds for any norm. For the second inequality, 
note that the matrix X must be positive definite hermitian and, in particular, non- 
singular and ||X-^|| = 1/A™„(X). By 



Hence 

iix-Mi = ^A^< 



CaU A, = \m^n{Xi) > 0. Then 

The right-hand side inequality holds true: A + A^^ > 2 for all A > and therefore 



Proof of Proposition: We have, as in psp . 



Pll<ll^n 11 ■ max ||X, -Xoll- ||X„||+ max || J^, || + 2||Xo | 
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From the previous Lemma we know that \\Xi\\ < M and \\X^ ^|| < for all < i < n. 
From Proposition |4] we have 

max \\Xi - Xi II = diam (a[Xo, Xj) < ( ■ d. 

Therefore 

||1|| < rn diam(A[Xo, ■ ■ ■ , X-^]) (M + M + 2M) < ( ■ 'idMrn? . 

\n + l J 

6 The infinite series for n = 3 

In this section, we expand Hamida's integral for n = 3 into a explicit infinite series 
with no integral operator, which can be directly implemented as a computer algorithm. 
Recall equation 

Tr / ((1 - tA)-^{dtA + tdA)f = 

J M 

Ja^ mi + ma + 3 ^ ^ 

mi >m2 >U 

Write A = Uo + siUi + S2U2 as in ([9]). Fix mi > ma > 0. We have 
Tr / 

= Tr / ((7o + siUi + s2U2y"^dA{Uo + s^Ui + S2U2T^dA 
Ja2 



= Tr 



/ 51/ 'Si^Sa^^ui%2 



where Bp ^ (f" > P + ?) is the sum of all matrices of the form Ui-^ . . . Ui^, ij £ {0, 1, 2}, 
containing p times f/i and q times U2 (and hence r — p — q times Uo). Reordering and 
using that dA — dsiUi + dsaC/a we get 



U<~Ui +li2 ^ "I'l 

E sr+"^5r+''^B™%2t^2S™^„,C/i dsidsa J . 



0<iJi+i;2^^2 

The second term can be rewritten, by a cyclic permutation, to have the form 
B^^%^UiB^^\^^U2 dsids2. That is, if we write 

-Em,,! = E Si^+"^S2^'''"2B^j_„2C/ii?"j_„2i7a dsidsa 

0<lli+U2<^^ 
0<l'l +1^2 

then (|17p equals Tr /_^2 (-E'mi,m,2 — ^'m2,™i)i for each mi > ma > 0. Hence for each pair 
of arbitrary integers mi, ma > we have Emi,m2 when mi > ma and —Emi,m2 when 
mi < ma. Consequently, if we substitute (mi — ma) by |mi — ma|, we can write p6|l 
as 

3Tr/ Y: (18) 
Ja2 mi + ma +3 

mi ,m2 >U 

The trace and the integral operator are hnear, and 
Tr f Em„m,= V Tr(BriU^^i^^%2t^2) / s'l'-^^' s^^^^^ ds,ds2 . (19) 

U<.Til+U2 5:^l 
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The integral can be calculated using integration by parts. If n, m > 0, 

/ S1S2 dsias2 ~ Si [ S2 ds2 fitsi — / — ^ dsi 

Ja-^ Jo \Jq J Jo m + 1 

-dSl — . . . = / / ' —^dSl 

Jo 



(m + l)(m + 2) Jo (m + l)...(m + n + l) 



(m + 1) . . . (m + n + 2) (7i + m + 2)!' 

The calculation is still valid when n — and, by symmetry, when m = 0. 

On the other hand, each term B'!^^^^^U-iB^^„^U2 amounts to the sum of matrices 
in the form (7i^ . . . Ui^_^ UiUj^ . . . Uj^^ U2 with ui respectively U2 matrices Ui resp. U2 
among the first mi, and vi resp. V2 matrices Ui resp. U2 among the last m2- By cyclic 
permutation, we can write it as 

UlUj,...Uj^^U2U^,...Ui^^ . (20) 

Given A; > 1 and an array UiUi^ . . . Ui^, how many times (if any) does it appear 
in (|18p and with which coefficient? Compare the expression UiUi^ . . . Ui^. to H20p . For 
every j with ij — 2 we have a product as in (|20p with 



m2 


= J-1 






nil 


= k-j 






Vl 


= \{ip = 


f 


11<P< j-l}| 


V2 


= \{ip = 


2 


|1<P< j-l}| 


Ul 


= \{ip = 


f 


\j + l<p<k}\ 


U2 


= \{ip = 


2 


j + 1 <p< fc}| 



Define ni = ui + vi and 712 = U2 + V2, that is, 

m = l{jp = l|l<p<A:}l 

n2 = l{jp = 2 1 f < p < fc}| - I . 

Let us denote by X2 the characteristic function on 2 (X2(2) = I, X2(*) = if i 7^ 2). 
Using (|19p and the previous calculations, we can rewrite (|18p as 

00 2 

fc — 1 

where 

E*" /. w ,xfe|fc-2j + f| (iti + ui)!(tt2 + U2)! 
. fc + 2 (ni+.i+.2+^.2 + 2)! 

(-f)'' (ni)!(n2)! ^ , 
= m (m+n2 + 2)! Ex.fe)|fc-2j + l|- 

Note that we have use that, for positive definite hermitian matrices Xi of size n and 
the spectral norm, Tr (J^^^^ X„) = Er=i Tr(X„), since < Tr(X) and Tr(X) < 
n\\X\\. 

Proposition 14 (Infinite series for n — 3) // ||A|| < I then 
Tt [ {{l-tAy\dtA + tdA)f = 

J M 

fc=i n,...,ifc=o 
where ni, n2 and c depend on each tuple (ii, . . . ,4^4.1) as 

ni = \{ip = 1 1 1 < p < k}\, 
na = |{jp = 2|l<p<fc}l-l, 

c = ^X2(*.)|fc-2j + l|. 
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7 Computacional aspects and algorithm 



We discuss the computational aspects of our algorithm. We write || || for the spectral 
norm. 

7.1 Initial steps 

The program is given matrices Yo,Yi,Y2,Y3 £ GLn{C) (normally N — 7). It first 
computes the positive definite hermitian matrices Xq, Xi, X2, such that X*Xi = 
Y*Yi for all < j < 3 (see Remark [2]|. This is done via the Schur decomposition of 
X* Xi. Then the following values are stored 



dij — 1 1 Xi — Xj 1 1 — \/ \max {{Xi — Xj)*[Xi — Xj)) , 
f^i = ll^ill ~ ^max{Xi)^ 

= ii^r'n = i/A™.„(xo, 

for i,j G {0,1,2,3}. 

7.2 Barycentric subdivision 

This is a recursive step. Given a 4-simplex A[Xo, . . . , X3] of the fc"'-barycentric sub- 
division of A = A[Xo, . . . ,^3] (just the initial simplex A when fc = 0), we define the 
corresponding matrix A as in Proposition 1121 and try to prove that ||A|| < 1 using 

Pll < ll^u'f • max -Xoll- max || + 2||Xo | 

If we succeed, we move to the next step. If we fail, we calculate the barycentric 
subdivision of A[Xo, . . . , X3] and apply this step recursively to each of the resulting 24 
simplices. 

Remark 15 Provosition M'A guarantees that the recursion finishes: the maximal depth 
required in the barycentric subdivision is bounded above by j log^^^ij^fi^gl^) | • 

7.3 Infinite Series 

We are given (Xq, . . . ,Xn) with the associated matrix A satisfying ||j4|| < 1. Our aim 
is to approximate the value of the infinite series in Proposition [TJ] We keep the same 
notation as there. Also, let us write fii...i^. for the factorial coefficient, that is, 

, _ (m)!(n2)! 



(m +n2 +2)!' 

The [k + l)th term in the series (|21|l can be obtained from the fcth term as follows. For 
each ii, . . . ,ik G {0, 1, 2}, we have three new matrices 

{U^U,,...U,,)U, j = 0,l,2 

and three new factorial coefficients 

fii - iki = \ f. ^ ^ "j _ 1 9 (22) 

where ni and n2 are the ones defined for the tuple (ii, . . . , ifc). Hence for each such 
tuple (ii, . . . ,ik), we should store the values of UiUi^ . . . Ui^ and fii...ik- The other 
terms in the series can be inexpensively calculated in each iteration, see (|21|l . 

7.4 Error estimate 

It remains the question of bounding the error after each iteration. We now verify the 
absolute convergence of the series (|21|l with respect to the spectral norm by finding a 
majorising series, which serves as means of estimating the error after each iteration. 
If X is a matrix with eigenvalues Ai, . . . , A„ then 



|Tr(X)l = ^A, <n\Xrr,ax\. 
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So if X is hermitian, \Tr{X)\ < n\\X\\. Write A ^ Uo + siUi + S2U2 as in Sectiong] 
Suppose that \\Ui\\ < m, i — 0,1,2 (we will need m < 1/2). Then 

|Tr ([/iC/,, ...U^,)\<4 \\UiU,, . . . (7,J| < 4m'=+' . 

On the other hand 



E |fc-2j + ll = { 



/2 if fc even, 

(fc2-l)/2 iffcodd, 



so we have, 

Additionally, one can prove that 



(?ii)!(tt-2)! ^ 2^-1 

I'm -I- 4- 9^ ~ 



il,...,ifc=0 

by induction on k: for A; = 1 we get 5/6 < 1 and the induction step follows from 
which implies 

2 2 

fil---ik + l — 5^ 2/ii...ij. < 2 . 

All in all, if we call the fcth coefficient of the series (I21II , we have 



2 



l«*l = fcT2 ^ /n.....y|Tr(f/iK,...[/, 

il,...,ifc=0 



< 



<^(2mr^ 

Therefore, the series (|2ip is majorised by the convergent series 



if m < 1/2. Hence we need each \\Ui\\ < m < 1/2, as expected, since this immediately 
implies < 1. 

Let Sn = Ofc be the value of the series (|2ip after n iterations and S = X^fcLi "^fc 

the exact value. Since (I23|l is a majorising series, the error after n iterations \S — Sn\ 
is bounded above by 

^ 2 ^ ^ 2 ^ ^ (1 - 2m)2 ^ 2 ^ ^ 

fe=i fc=i ^ ' fc=i 

8 Test cases 

We have implemented an algorithm in Matlab 13^ to compute the series in Proposition 
1141 We have performed some test cases where we compare the program's output with 
the expected value of Hamida's function tp. 
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8.1 Test case 1 



Let d be a prime, u £ Q(\/— d)* and n,m £ Z. Set 

f u \ f 1 m \ . /I ny/^ \ 

3'=[o u )' ''-[o I ) '''"^ 3'-[o 1 )■ 

Let 

z = (1,51,5152,515253) - {i, 91,9193,919293) + (1,53,5153,515253) 
-(1, 52, 5152, 515253) + (1, 52, 5253, 515253) - (1, 53, 9293,919293)- 

Then z is a 3-cycle and ifiiz) — \Sl. Write the cycle as z = ci — C2 + C3 — C4 + cs — ce • We 
obtain answers close to zero for Hamida's function already at each Ci, that is, ifiici) ~ 
for each 1 < i < 6. The reason is that Hamida's function is unitarily normalized. If 

U2 = 1 J and «3 = ^ ), 

then M252 and M353 are the identity, and every Ci above can me made into a cycle 
with two consecutive repeated elements. In this case v'(ci) ~ 'p{Ci) — 0, since Hamida's 
integral vanishes for any chain with two repeated values (an n-form in n — 1 variables 
is zero). 

The algorithm fails to detect the vanishing directly but gives answers close to zero. 
We keep track of the positive respectively negative partial sums (the sum of the positive 
respectively negative terms of the series in Proposition I14p so that the answer is close 
to zero up to an error relative to the order of the positive and negative sums. 

Additionally, we had to limit the maximum number of iterations to obtain an answer 
in a sensible time. Sometimes we modified the bound test (see Section [7. 2|) to force a 
deeper barycentric subdivision in order to compensate for the lack of iterations. 



Case 6? = 3, u — 1, m = l and n = 1 

In this case there are two repeated entries in each chain and we obtain answers close 
to the machine epsilon in our computer. 





POS 


NEG 


POS+NEG 


Cl 


79.212543798236482 


-79.212543798232005 


4.476 


10" 


12 


C2 


35.255902765497538 


-35.255902765496728 


0.810 


10" 


12 


C3 


25.761457824702696 


-25.761457824702926 


-0.231 


10" 


12 


C4 


22.698322652159643 


-22.698322652159863 


-0.220 


10" 


12 


C5 


35.423517445833234 


-35.423517445834165 


-0.931 


10" 


12 


C6 


38.702788638164705 


-38.702788638165714 


-1.009 


10" 


12 



Case rf = 3, u = 2, TO = 1 and n = 1 

This time the algorithm was not simplified by repeated entries. We had to forced an 
increased barycenter subdivision and limit the iterations to 4 to finish in a reasonable 
time. The output is summarized in the following table. 





POS 


NEG 


POS+NEG 


Cl 


402.6517675 


-402.6512473 


0.0005202 


C2 


778.9142690 


-778.9130614 


0.0012076 


C3 


933.8326019 


-933.8463496 


-0.0013748 


C4 


489.5582685 


-489.5751380 


-0.0168695 


C5 


494.6188740 


-494.6213859 


-0.0025119 


C6 


449.8738903 


-449.8708292 


0.0030611 



The solutions are nevertheless close to zero relative to the positive and negative values. 
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Case d = 3, u = a/— 3, to = — 1 and n = 2 

The same remarks apply to this case. The output is given again in table format. 





POS 


NEG 


POS+NEG 


Cl 


1390.3704696 


-1390.3769837 


-0.0065141 


C2 


4185.0669067 


-4185.0619767 


0.0049300 


C3 


5879.8887455 


-5879.9673981 


-0.0786525 


C4 


1701.0187054 


-1701.0253901 


-0.0066847 


C5 


2774.6657137 


-2774.6719721 


-0.0062584 


C6 


2601.1767976 


-2601.1739894 


0.0028082 



As before, the answers are close to zero with small error relative to the positive and 
negative values. 

8.2 Test case 2 

Let F be a field and u£ F with u" = 1. Set 

/mOO\ ( u Q Q \ I 1 0\ 

o= 0u-^0 ,fe= 1 ,m= v\. 
\0 1/ \00m-^/ \0 -y-^ / 

Note that ah = ba, a" = e, waw^^ — b and wbw~^ = a. Define (in bar resolution 
notation) 

zi = [w\a\b] - [a\b\w] - [w|6|a] + [b\a\w] + [a\w\a] - [b\w\b], 

n-l 

Z2 = [e|e|&] + [b|e|e] + y ^ [a \a\b] — [a''|6|a] + [b\a \a]. 

r = l 

Then d{zi) = 2{a,b} and d{z2) = n{a,b} so [nzi - 2z2] G H3 {SL3{F);Z) [5]. If F is a 
number field, this class should be in the kernel of the Borel regulator. 

We have done several test cases for F = Q{u), u = g^'^'''" and v — 1. In all cases, 
no barycentric subdivision is needed (the matrices are already close enough) and the 
algorithm takes just a few seconds to finish (hence no need in general to limit the 
number of iterations). Moreover, the values of ip at each chain forming z\ and Z2 are 
deeply related, as we show. 

To facilitate reference we write = ci — C2 — C3 + C4 + cs — ce and 22 = di + d2 + 

da — d4 + rfs + • ■ • + rf3(n-l)+2- 

Case n = 3 

The results are summarized in the following two tables. For once, we show the positive 
and negative answers, altough this is not needed anymore. 





POS 






NEG 






POS-fNEG 




Cl 





















C2 


0.54738221262 


10" 


-47 


-0.18246073754 


10" 


-47 


0.36492147508 ■ 10" 


-47 


C3 





















C4 


0.54738221263 


10" 


-47 


-0.18246073754 


10" 


-47 


0.36492147508 • 10" 


-47 


C5 


0.54738221263 


10" 


-47 


-0.18246073754 


10" 


-47 


0.36492147508 • 10" 


-47 


C6 


0.54738221263 


10" 


'47 


-0.18246073754 


10" 


-47 


0.36492147508 ■ 10" 


-47 


Zl 
























POS 






NEG 






POS-I-NEG 




di 























d2 


0.54738221262 


10" 


-47 


-0.18246073754 


10" 


-47 


0.36492147508 


10" 


-47 


dz 


0.54738221262 


10" 


-47 


-0.18246073754 


10" 


-47 


0.36492147508 


10" 


-47 


di 


0.54738221262 


10" 


-47 


-0.18246073754 


10" 


-47 


0.36492147508 


10" 


-47 


ds 


0.54738221262 


10" 


-47 


-0.18246073754 


10" 


-47 


0.36492147508 


10" 


-47 




0.32842932758 


10" 


-46 


-0.10947644252 


10" 


-46 


0.21895288505 


10" 


-46 


d7 


0.27369110631 


10" 


-46 


-0.07298429502 


10" 


-46 


0.20070681130 


10" 


-46 


ds 


0.21895288505 


10" 


-46 


-0.05473822126 


10" 


-46 


0.16421466379 


10" 


-46 


Z2 


0.38316754884 


10" 


-46 


-0.12772251628 


10" 


-46 


0.25544503256 


10" 


-46 
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Note that, if x = 0.364921475084588 ■ 10""^, then da = 6x, dr = ^-x and ds = fa;. 
Write c = {(p{ci), . . . , 'p{ce)) and d = {(p{di), . . . , (^(ds))- Then we have showed that for 
n = 3, c = (0,x, 0, x,x, x) and d = (0, x,x, x, x, 6x, ll/2a;,9/2a;). 



Case n = 4 

We obtain ci = . . . = ce = di = . . . = dn = all zeroes, including the positive and 
negative subtotals. 



Case n = 5 

The output for n = 5 is 









Cl 









d5 


0.36492147508 


10" 


-47 


C2 


0.36492147508 


10" 


-47 


d(i 


-0.91230368771 


10" 


-47 


c-i 









d7 


-0.91230368771 


10" 


-47 


C4 


0.36492147508 


10" 


-47 


ds 


-0.18246073754 


10" 


-47 


C5 


0.36492147508 


10" 


-47 


do 


-4.01413622593 


10" 


-47 


C6 


0.36492147508 


10" 


-47 


dio 


-0.72984295017 


10" 


-47 


di 









dii 









d2 


0.36492147508 


10" 


-47 


di2 


-5.47382212627 


10" 


-47 


d3 


0.36492147508 


10" 


-47 


di3 


-5.47382212627 


10" 


-47 


d4 


0.36492147508 


10" 


-47 


di4 


0.36492147508 


10" 


-47 



That is, for n = 5 we obtain 

c = {0,x,0,x,x,x) and 

d = (0, K, x,a;,a;, — 5/2a;, — 5/2a;, — l/2a;, — 11a;, — 2a;, 0, — 15a;, — 15a;, a). 



Cases n < 11 

We always obtain results which are rational multiples of x. We summarize the infor- 
mation in the following tables. 



n 




c 




3 




{0,x,0,x,x,x) 





4 




(0,0,0,0,0,0) 





5 




(0, X, 0, x,x, x) 





6 




(0, a:, 0, a:, 0, 0) 





7 




(0,0,0,0,0,0) 





8 




(0,0,0,0,0,0) 





9 


(0, 


-l/2a;,0, -l/2a;,a;,a;) 





10 




(0,a;,0,a;,0,0) 





11 




(0, a;, 0, a;, 0, 0) 






d ^{^2) 

(0,0,0,0,0,0,0,0,0,0,0) 

cc ^ cc ^ «Z/j X ^ ^ ^ 2' J 1 111/ ^ 2^ ^ 0^ 1 5tj/ ^ 153/ ^ x"^ 2~ ^ 

(0, X, fa;, a;, x, 2x, 2x, x, ^a ^T-a; 

(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -4a;, -4a;, 0, -10a;, -4a;, 0, -24a;, -20a;, 0) -lOx 

(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) 

-191a;, 11a;, -205a;, -354a;, 13a;, -370a;, Sf-x, fa;, -622a;, -598a;, 13a;) ^x 

^0 ^ X ^ X ^ X ^ X ^ X ^ X ^ X ^ ^ X ^ X ^ X ^0^0^ X ^ ^ x ^ ^ x ^ ^ x ^ 

(0, x,^x,x,x,^-x, 8x,x, 15x, fx, x, 16x, 16x, 1, 48x, 40x, x, 54a;, 

108a;, a;, 231a;, 108a;, x, 126a;, 245a;, x, 234a;, 224a;, a;, 243a;, 243a;, x) - fa; 

The fact that all the values obtained are rational multiples of one another is very 
significant. We have predicted that rvp{zi) — 2ip{z2) = but instead we obtain that 
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ifiizi) always vanishes and ip{z2) 7^ 0. We suspect that we have made a minor mistake 
at some stage, since we are in each case very close to get ip{z\) = 2x and ifi{z2) — nx, 
with a readjustment of the signs. 
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